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It is shown that different pairs of stress-energy and spin tensors of quantum relativistic fields 
related by a pseudo-gauge transformation, i.e. differing by a divergence, imply different mean 
values of physical quantities in thermodynamical nonequilibrium situations. Most notably, transport 
coefficients and the total entropy production rate are affected by the choice of the spin tensor of the 
relativistic quantum field theory under consideration. 



I. INTRODUCTION 

In a previous paper [l[ we showed that in relativistic quantum field theory, different pairs of stress-energy and spin 



C/3 



> 



tensors related by a pseudo-gauge transformation 0, 0| 



(N ; 

where $ is a rank three tensor field antisymmetric in the last two indices (often called and henceforth referred to 
as superpotential) and Z a rank four tensor antisymmetric in the pairs aX and fj,v, are indeed thermodynamically 
incquivalent as they imply different mean values of physical quantities for a rotating system at equilibrium. Particu- 
^■i larly, for the free Dirac field, we showed that the canonical and Belinfante (obtained from the canonical one by setting 
$ = S and Z = in (^Q), hence with a vanishing new spin tensor S') quantum stress-energy tensors result in different 
mean values for the momentum density and the total angular momentum density. 

The thermodynamical inequivalence is (at least in our view) surprising because it was commonly believed that the 
only physical phenomenon which can discriminate between stress-energy tensors of a fundamental quantum field theory 
^v^j ' related by a transformation like ([1]), preserving total energy, momentum and angular momentum, is gravity, or, in other 
words, the coupling to a metric tensor. In this paper we reinforce our previous finding by showing that the inequivalence 
{SJ ■ extends to nonequilibrium thermodynamical quantities, specifically entropy production and transport coefficients. 
^sO \ In summary, we will show that the use of different stress-energy tensors, related by (JlJ, to calculate transport 
■ coefficients with the relativistic Kubo formula leads, in general, to different results. Therefore, at least in principle, 
an extremely accurate measurement of transport coefficients or total entropy in an experiment where dissipation 
is involved, would allow to disprove a candidate stress-energy or spin tensor, with obvious important consequences 
in relativistic gravitational theories. This finding means, in other words, that the existence of a fundamental spin 
tensor affects the microscopic number of degrees of freedom, or at least on how quickly macroscopic information gets 
• *h ■ converted into microscopic, namely on entropy generation. 

The paper is organized as follows: in Sect. [TT] we will extend the framework of the nonequilibrium density operator 
introduced by Zubarev [ij to the case of a non- vanishing spin tensor. In Sect. lIIIl it will be shown that the nonequilib- 
rium density operator is not invariant under a pseudo-gauge transformation fl} , that is it does depend on the chosen 
couple of stress-energy and spin tensor. In Sect. ITVl we will provide a general formula for the change of mean values of 
observables and we will determine how entropy is affected by a pseudo-gauge transformation. In Sect. |V] we will show 
that transport coefficients are also modified and, particularly, we will focus on the modification of the Kubo formula 
for shear viscosity. Finally, in Sect. IVI1 we will discuss the implications of this finding and draw our conclusions. 



Notation 



In this paper we adopt the natural units, with K= c = K = 1. 
The Minkowskian metric tensor is diag(l, —1, —1, —1); for the Levi-Civita symbol we use the convention £ 0123 = 1. 
We will use the relativistic notation with repeated indices assumed to be saturated. Operators in Hilbert space will 
be denoted by an upper hat, e.g. R, with the exception of the Dirac field operator which is denoted with a capital 
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II. NONEQUILIBRIUM DENSITY OPERATOR 

A suitable formalism to calculate transport coefficients for relativistic quantum fields without going through kinetic 
theory was developed by Zubarev 0, 0] , extending to the relativistic domain a formalism already introduced by Kubo 
Q. In this approach, a non-equilibrium density operator is introduced which reads Q 1 : 



p = iexp[-T] = -^cxp 



/■*' 

lim e 



dte £ (*-^ f d 3 x (f 0v (3 v {x)-3°£,{x)) 



(2) 



where j is a conserved current, the four-vector field /3 is a point-dependent inverse temperature four-vector (/3 = u/Tq, 
u being a four- velocity field and To the comoving or invariant temperature) and £ = Po/Tq a scalar function whose 
physical meaning is that of a point-dependent ratio between comoving chemical potential po and comoving temperature 
To; the Z factor is analogous to a partition function, i.e. a normalization factor to have trp = f . The operators in the 
exponential of Eq. ([2]) are in the Heisenberg representation. It should be stressed that in the formula ([2]) covariance 
is broken from the very beginning by the choice of a specific inertial frame and its time. However, it can be shown 
that the operator p is in fact time- independent Q, namely independent of t', so that p is a good density operator in 
the Heisenberg representation. 

In the formula ([2]) the possible contribution of a spin tensor is simply disregarded; therefore, the formula is correct 
only if the stress-energy tensor is the symmetrized Bclinfante one (or improved ones, see last section), whose associated 
spin tensor is vanishing. It is the aim of this Section to find the appropriate extension of the formula ^ with a spin 
tensor. 

Using the identity: 

e e(t-t') (fOu Mx) = (±. e ^l\ (T^P V {X) -?£(*)) 

integrating by parts and taking into account the continuity equations d^T^ v — — 0, the operator T in Eq. @ 
can be rewritten as: 



T = 



(T 0l %(t',x) -?£(f,x)) + lim f dt e 6 '*-*') f dSn, (f iv p v {x) - j*f(s)) 



- lim f* die 6 **"*') j d 3 x (T^dMx) -i^tfa:)) (3) 



The first term the so-called local thermodynamical equilibrium one, which is defined by the same formula of the global 
equilibrium [|| with x-dependent four-temperature and chemical potentials, whereas the term dependent on their 
derivatives is interpreted as a perturbation. 

At equilibrium, the right hand side should reduce to the known form, which, at least for the most familiar form of 
thermodynamical equilibrium with f3 cq = (I/T, 0) = const and £ oq = p/T = const is readily recognized in the first 
term setting (3 = f3 cq and £ = £ oq : 



d 3 x (f °X q - j°f q ) + lim j dt e 6 ^*') j dS m (f lu (3^ - ?£ cq ) 
- limy"* die^"*') j d 3 x (f^d^-fd^) =H/T-pQ/T 
+ lim / die'**-*') f dSn, (f w p cq - j ,; £ oq ) - lim f dt e B <*-*'> j d 3 x (f^d^ cq -J^" 5 ) (4) 



Hence, the two rightmost terms of Q must vanish at equilibrium. Indeed, the surface term is supposed to vanish 
through a suitable choice of the field boundary conditions while the third term vanishes in view of the constancy of 



Throughout the paper, the four-vector x implies the time t and position vector x, i.e. x = (t, x). The dependence of the stress-energy 
and spin tensor on x will always be understood. 
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f3 cq and £ cq . However, this is not the case for the most general form of equilibrium; in the most general form (see 
discussion in ref. [6|), whilst the scalar £ eq stays constant the four- vector j3 fulfills a Killing equation, whose solution 
is 0: 

/3?(x) = b? + U %x» (5) 
with both the four-vector b eq and the antisymmetric tensor w oq constant. Therefore: 

which in general is non- vanishing, so that the third term on the right hand side of Eq. (|4]) survives. For instance, for 
the thermodynamical equilibrium with rotation [6|, the tensor ui turns out to be: 

w£ = «/r (6) 

u> being the angular velocity and T the temperature measured by the inertial frame. 

In order to find the appropriate generalization of the operator T, let us plug the formula ([5]) of general thermody- 
namical equilibrium into the (@J: 



x (f 0u /3^ - j°£ cq )) + lim I dt e £ (*-*') j dSn, (f lv {b cq + w° q ^) -?£ cq ) 



+ lim J dt e^*"*') J d 3 x f ""wJJ (7) 

where d /J £ cq = has been taken into account. For a symmetric stress-energy tensor T, the last term vanishes, but 
if a spin tensor is present T may have an antisymmetric part. Particularly, from the angular momentum continuity 
equation: 

T^u* = \ (T^ - = -\d x S x ^u° q (8) 

so that the last term on the right hand side of Eq. (Jjj can be rewritten as: 

i-t' e 1 ft' 



lim f die^*"*' 5 [ d 3 xf u V5 = -~w!£lim / dt e 5 ^^ f d 3 x dxS^ 

io;° q lim y d 3 x J dt e^-^^S '^ - ^uj cq , lim J dt e<*~^ J dSn.S 1 '^ 



lim / (]- : x / 

The first term on the right hand side of © can be integrated by parts, yielding 



1 -ojH Vim J d 3 x £ dt e^-O^o,^ = J d ^s^^t\ x ) + ^l \hne £ dt e'W J d 3 xS°^ 



(9) 



(x) (10) 



Plugging the Eq. (fit))) into ([5]) and this in turn into ([7]) we obtain 



1 TH3 



T oq = y d 3 x (t^^-j ^-^^'^ +lim y dte £ (*-*') b cq J dSn t T w -£ cq J dSS 

- -lo2 [ dSni{x^f iv -x v f" i +§ 1 '^)] +iw° q lime / dt e^ 4 "*') / d 3 x S°'" u (x) (11) 
2 F J J 2 ^ e-^o y 

where the surface term involving T in Eq. ([7]) has been rearranged taking advantage of the antisymmetry of the uj 
tensor. The surface terms in the above equations now are manifestly the total momentum flux, the charge flux and the 
total angular momentum flux through the boundary. All of these terms are supposed to vanish at thermodynamical 
equilibrium through suitable conditions enforced on the field operators at the boundary, so that the (jTTJ) reduces to: 



T oq = y d 3 x (f "^ - ?C cq ~ l^S '^ + l -uj cq lime J dt e^*"*' 5 J d 3 x5 D ^(a 



(12) 
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The first term on the right hand side just gives rise to the desired form of the equilibrium operator. For instance, for 
a rotating system with uj as in Eq. ([5]) one has [6| : 



d 3 x (f ov ^-Te q -l^s°^A 



H/T — fxQ/T — ujJ/T 



J being the total angular momentum, which is the known form [TT(. Nevertheless, the second term in Eq. (fT2|) does 
not vanish and, thus, must be subtracted away with a suitable modification of the definition of the T operator. The 
form of the unwanted term demands the following modification of @ : 



~ 1 r 1 

P = 7?cxp[-TJ = -exp 



— lim s 

£->0 



dt J d 3 x (f 0v p v {x) - is^uvOr) 



(13) 



where uj^ v (x) is an antisymmetric tensor field which must reduce to the constant uJ;J tensor at equilibrium. It is easy 
to check, by tracing the previous calculations, that the equilibrium form of T reduces to the desired form: 



T cq = J d 3 x ffOi/geq _-j0^n _ ^•q^O./u^ 



as the spin tensor term in Eq. (JT2J) cancels out. Therefore, the operator (ITS")) is the only possible extension of the 
nonequilibrium density operator with a spin tensor. 

The new operator T can be worked out the same way as we have done when obtaining Eq. ([3]) from Eq. ([2]): 

f = [ d 3 x (f Q ^ v {t'^)-3 Q i{t'^)-\^ v ^ v {t'^ 



+ lim J* dt e 6 **-*') J dSn t (f iv /3 v (x) -?£(x) ~ -^'""^(a:)) 



\ lim J dt e £ (*-*') f d 3 x (f^(d^(x) + d^ v (x)) + f^(d^(x) - d^ v {x) + 2u MV (x)) 



(14) 



where: 



and the continuity equation for angular momentum has been used. The first term on the right hand side is the new 
local thermodynamical term whilst the third term can be further expanded to derive the relativistic Kubo formula of 
transport coefficients (see Appendix A). 



III. NONEQUILIBRIUM DENSITY OPERATOR AND PSEUDO-GAUGE TRANSFORMATIONS 

A natural requirement for the density operator (1131) would be its independence of the particular couple of stress- 
energy and spin tensor, because one would like the mean value of any observable O : 

O = tr(pO) 

to be an objective one. In ref. [l| we showed that even at thermodynamical equilibrium with rotation this is not the 
case for the components of the stress-energy and spin tensor themselves because they change through the pseudo-gauge 
transformation ([T]). However, at equilibrium, p itself is a function of just integral quantities (total energy, angular 
momentum, charge) which are invariant under a transformation ([1]) provided that boundary fluxes vanish, so a specific 
operator O, including the components of a specific stress-energy tensor, does not change under ([T]). However, it is not 
obvious that this feature persists in a nonequilibrium case, in fact we are going to show that, in general, this is not 
the case. 
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Let us consider the operator T in (fl3|) and how it gets changed under a pseudo-gauge transformation ([TJ with 
Z = 0. The new operator T' reads: 



T' = X + ~ lim e y" di e^"*') J d 3 x (dxft '" p u (x) + ^'""ovO*)) 



where: 
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(15) 



(16) 



is antisymmetric in the first two indices. We can rewrite Eq. ([T5 



T' - T = - lim £ / dt / d 3 x e^*"*' 5 



as: 



-oo •/ 

1 /■*' 

- lim e / dt e £(t_ * 



dSm ft°' v p„{x) - J d 3 x {^'"dx^u - S '^^^)) 



after integration by parts. Let us now write the general fields /3 and u> as the sum of the equilibrium values 
perturbation, that is - 



P(x) = /3 cq (x) + 5/3(x) w{x) = w eq + Su(x) 



(17) 
and a 

(18) 



and work out first the equilibrium part of the right hand side of Eq. (|17p . As d\(3° q — —^° Xv one has: 



(f'-T)| eq = ilims / dte^ 



= I lim e / dt e £ (*-*') 



— OO 



dSn, (p l ^p cq {x) + J d 3 x ($ X ' 0l, u>H - S '* 1 ^ - &> xo ull + P'^iv^,) 



= — lim 



im£ / dte 5 ^-^ / dS ni<^°' u Pl^ix) 



(19) 



where we have used the Eq. (fT5|) and the antisymmetry of indices of the superpotential $. By using the Eq. the 
last expression can be rewritten as: 



fee / die^*"*' 5 



The two surface integrals above are the additional four-momentum and the additional total angular momentum, in 
the operator sense, after having made a pseudo-gauge tranformation ([T]) of the stress-energy and spin tensor. If the 
boundary conditions ensure that the momentum and total angular momentum fluxes vanish (in order to have conserved 
energy and momentum operators) for any couple (T, S) of tensors, then the two fluxes in the above equations must 
vanish as well. Therefore, we can conclude that: 



T Icq — T| cq 



Now, let us focus on the nonequilibrium perturbation of the T operator. 



(T'-T)| non _ oq = -lime / die £ (*-*') 



dSm ft°' u 5p v - ( d 3 x xo ' v d x 60 u - S^Sunv 



= - lim £ / dt e^*-^ 

2 e^O ' 



— oo 
4' 



= - lim £ dt e E{t ~ l,) 

2 e^O ' 



dSn, ^°' v 5p v - f d 3 x ($ A < 0,y - & M - $ v ' X0 )d x 6p v - P'^Su:^ 



dSn, p°-»Sf3„ " /d 3 x $ x >°»(d x 60 u + dj(3 x ) - S ^ Q(0 x J& - dj(3 x ) + 6lj Xv 



(20) 
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where the dependence of 5/3 and 8ui on x is now understood. It can be seen that it is impossible to make this difference 
vanishing in general. One can get rid of the surface term by choosing a perturbation which vanishes at the boundary 
and the last term by locking the perturbation of the tensor oj to that of the inverse temperature four- vector: 

Su> Xv {x) = ~(d x Sp v (x) - d v S/3 x (x)) (21) 

but it is impossible to cancel out the term: 

ST = I lim e [ it e £ (*-*') [ d 3 x &> 0u (d x 5/3Jx) + d u 5/3 x (x)) (22) 
2s->o 7_ oc J 

unless in special cases, e.g. when the tensor $ is also antisymmetric in the first two indices. 

We have thus come to the conclusion that the nonequilibrium density operator does depend, in general, on the 
particular choice of stress-energy and spin tensor of the quantum field theory under consideration. Therefeore, the 
mean value of any observable in a non-equilibrium situation shall depend on that choice. It is worth stressing that this 
is a much deeper dependence on the stress-energy and spin tensor than what we showed in ref. [l| for thermodynamical 
equilibrium with rotation. Therein, mean values of the angular momentum densities and momentum densities were 
found to be dependent on the pseudo-gauge transformation ([1]) because the relevant quantum operators could be 
varied, but not because the density operator p was dependent thereupon. In fact, at non-equilibrium, even p varies 
under a transformation ([1]). Note that, in principle, even the mean values of the total energy and momentum could 
be dependent on the quantum stress-energy tensor choice although boundary conditions ensure, as we have assumed, 
that the total energy and momentum operators are invariant under a transformation (TTJ). Again, this comes about 
because the density operator is not invariant under ([1]), in formula: 

tr(pP'") = tttfPi*) + tr(pP^) 

It must be pointed out that the variation of the Zubarev non-equilibrium density operator (|22p depends on the 
gradients of the four-temperature field and it is thus a small one close to thermodynamical equilibrium. In the 
next Section we will show in more details how the mean values of observables change under a small change of the 
nonequilibrium density operator, or, in other words, when the system is close to thermodynamical equilibrium. 

IV. VARIATION OF MEAN VALUES AND LINEAR RESPONSE 

We will first study the general dependence of the mean value of an observable O on the spin tensor by denoting 
by ST the supposedly small variation, under a transformation ([1]), of the operator T. This can be either the one in 
Eq. (|22|) or the more general (only bulk terms) in Eq. ((20)) . We have: 

MP 6) = ±jtT(exp[-T'}d) = i 7 tr(cxp[-f - ST}0) (23) 

being Z' = tr(exp[— T — ST]). We can expand in <5T at the first order (Zassenhaus formula): 
Z' ~Z- tr(exp[-T]<5T) 

tr(exp[-T - ST}6) ~ tr (exp[-¥ }(I - ST + i[T, ST] - i[T, [f, ST}] + . . .)&) (24) 

hence, with ( ) = tr(p ), at the first order in ST: 

tr@6) = (O)' ~ (6>(1 + (ST)) - (OSf) + i([T, ST}0) - i([T, [T, ST}}0) + ... 

which makes manifest the dependence of the mean value on the choice of the superpotcntial $. 

As has been mentioned, close to thermodynamical equilibrium, the operator ST is "small" and one can write an 
expansion of the mean value of the observable O in the gradients of the four-temperature field, according to relativistic 
linear response theory Q. This method, just based on Zubarev's nonequilibrium density operator method, allows 
to calculate the variation between the actual mean value of an operator and its value at local thermodynamical 
equilibrium for small deviations from it. In fact, it can be seen from Eq. (|22p that the operator <5T, from the linear 
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response theory viewpoint, is an additional perturbation in the derivative of the four-temperature field and therefore 
the difference between actual mean values at first order turns out be (see Appendix A for reference): 

A(6) ~ lim -e I At e e (*-''> / d 3 x ([$ A <°», d}) (d x 6/3„(x) + d„6p x {x)) (25) 
where (. . .)o stands for the expectation value calculated with the equilibrium density operator, that is: 

po = ^cxp[-£/T + A 4/T] (26) 

Since tr(po[ ( I )A ' 0l/ , O]) = tr($ A ' 0l/ [O, po]) the right hand side of (|25|) vanishes for all quantities commutating with the 
equilibrium density operator, notably total energy, momentum and angular momentum. Nevertheless, in principle, 
even the mean values of the conserved quantities are affected by the choice of a specific quantum stress-energy tensor, 
though at the second order in the perturbation 8/3. 

We now set out to study the effect of the transformation ([1]) on the total entropy. In nonequilibrium situation, 
entropy is usually defined as [Tlj the quantity maximizing — tr(plogp) with the constraints of fixed mean conserved 
densities. The solution ple of this problem is the local thermodynamical equilibrium operator, namely: 

exp[-/d 3 x (f^p^x) -3°Z(x) - lS°^^ u (x))} 

PLK(t) = ^— — (27) 

tr(exp[- J d 3 x [T^p v {x) - ~ §S°-^av(a:)J]) 

which - as emphasized in the above equation - is cxplicitely dependent on time, unlike the Zubarev stationary 
nonequilibrium density operator (|I3[) ; of course the time dependence is crucial to make entropy 

S= -tr(p LE logp LE ) (28) 

increasing in nonequilibrium situation. In order to study the effect of the transformation (TTJ) on the entropy it is 
convenient to define: 



t le = J d 3 x (f^^-Tm-^s '^^) 



(29) 



for which it can be shown that, with calculations similar to those in the previous section, the variation induced by 
the transformation (flj is: 



(30) 



As has been mentioned, it is possible to get rid of the surface and the last term through a suitable choice of the 
perturbations, in fact not of the second term. 

Since ^Tle is a small term compared to Tle we can determine the variation of the entropy (|28[) with an expansion 
in 5Tle at first order. First, we observe that (see also Eq. (|24p): 

Zl E = tr (cxp[-T LE - 5T LE ]) ~ tr (exp[-T LE ](J - 5T LE )) = Z LE (1 - (6T LE ) f ) 
where ( )^ stands for the averaging with the original Tl E local equilibrium operator. Hence, the new entropy reads: 
S' = -^tr(cxp[-T LE -<5T LE ](T LE + 5T LE )) +logZ{ jE 

LE V ' 

' (l + (^T LE ) ¥ )tr(exp[-T LE -,5T LE ](T LE + ( 5T LE )) + log Z LE + log(l - (ST LE > ? ) (31) 



•^LE 

We can now further expand the exponentials as we have done in Eq. (|24[) . First: 

tr (exp[-T LE - <5T L e]Tle) ~ tr ^exp[-T LE ](7 - <5T L e + T^f? le, ST le] - ^[Tle, [Tle,^Tle]] + ■ • O^le 
= tr(cxp[-T LE ]T LE ) - tr(exp[-T LE ]<5T LE T LE ) = Z LE (T LE ) ? - Z LE (<5T LE Y LE ) ? (32) 
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where, in the second equality, we have taken advantage of commutativity and ciclicity of the trace. Then: 

Irii -A , 1 



tr ^exp[— Tle - <5Tle](5Tle) — tr ^exp[-T L E](/ - <5T L e + -^[Tle, <5T le ] - - [T LE , [T LE , £T LE ]] + ■ • .)SYie 
~ tr(exp[-T L E]<5T LE ) = Z LE (ST LE )f (33) 
keeping only first order terms. Thus, Eq. (|3 1[) can be rewritten as: 

S' ~ -^(l + ( ( 5T L E> ? )tr(exp[-T L E-<5T L E](T L E + 5T LE )) + log Z LE + log(l - (ST LE ) 

^LE - 7 

1 



' T ) 

(1 + (5Tle)t ) (^le(Tle)x - Z LE (ST LE T le ) x + Zle(<5T L e) x ) + logZ LE + log(l - (<ST LE ) X ) 



•^LE 

= (1 + {ST LE ) f ) (<T LE } X - <<5T LE T LE ) X + (Sf LE ) f ) + log Z LE + log(l - <<5T LE } X ) (34) 

Retaining only the first order terms in <5Tle, expanding the logarithm for (5Yle)le *C 1 and inserting the original 
expression of entropy: 

S' ~ S - (8T LE T L e)t + (5T LE )f (Tle) x (35) 

Therefore, the variation of the total entropy is, to the lowest order, proportional to the correlation between T and 
(5T, which is generally non- vanishing. 

We can expand the above correlation to gain further insight. For the 5Tle, let us keep only the second term of the 
right hand side of Eq. (jBTij) : 

<5T LE = - J / d 3 x $ x ' »{d x 5l3 u + d v 5p x ) (36) 



2 . 

By using the (p9)l and the (|36|) . the Eq. (f35|) can be rewritten as: 
S'(i)~S(t) + \J d 3 x /" dV (f^{x)T°^))^-^ v {x))^^ 

- i|d 3 x| d 3 x' (<$ A ^(x)?(a ; 0) x -($ A ' Ov (x)} ¥ (i o (x0) ? )c(^0(9A^(^ 

- J/d 3 x / d 3 x' (($^(x)^(x0) x -<$ A ' o ^x)) ¥ (^(x0) ? )^(^0(^^ 

(37) 

where cc and a;' have equal times. The above expression could be further simplified by e.g. approximating the local 
equilibrium mean ( ) x with the global equilibrium one ( )o, but this does not lead to further conceptual insight. The 
physical meaning of Eq. (|37p is that the entropy difference depends on the correlation between local operators in two 
different space points multiplied by a factor which is at most of the second order in the perturbation 8(3. This kind of 
expression resembles the product of transport coefficients expressed by a Kubo formula times the squared gradient of 
the perturbation field. Therefore, the difference between entropies suggest that the introduction of a superpotential 
may lead to a modification of the transport coefficients. We will show this in detail in the next Section. 

V. TRANSPORT COEFFICIENTS: SHEAR VISCOSITY AS AN EXAMPLE 

As has been mentioned, a remarkable consequence of the transformation J!} is a difference in the predicted values of 
transport coefficients calculated with the relativistic Kubo formula, which is obtained by working out the mean value 
of the stress-energy tensor itself with the linear response theory and the nonequilibrium density operator in Eq. ©. 
For this purpose, the derivation in ref. @ must be extended to the most general expression of the nonequilibrium 
density operator including a spin tensor, that is, Eq. (113[) : it can be found in Appendix A. 

The equation (|25p. yielding the difference of mean values of a general observable under a transformation ([T]), cannot 
be straightforwardly used to calculate the mean value of the stress-energy tensor setting O — T^ u (y) because T^ v (y) 
gets transformed itself. It is therefore more convenient to work out the general expression of the Kubo formula and 
study how it is modified by ([!]) thereafter. 
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We will take shear viscosity as an example, the transformation of other transport coefficients can be obtained with 
the same reasoning. Shear viscosity, in the Kubo formula, is related to the spacial components of the symmetric part 
of the stress-energy tensor. It is worth pointing out that, since a non- vanishing spin tensor can make the stress-energy 
tensor non-symmetric, there might be a new transport coefficient related to the antisymmetric part of the stress-energy 
tensor. 

For the symmetric part of the stress-energy tensor Tg V = (l/2)(T A " y + T^'), using the general formula of rela- 
tivistic linear response theory (Eq. 170)) of Appendix A), the difference 5Tg V (y) between actual mean value and local 
equilibrium value reads, at the lowest order in gradients: 



5T^(y) 



l im L f df 1 e£(t _ / d 3 x ( T" a (x),T^(y) ) Q d p S^(x) 



— lim — 

2 e->o i 



t' 



dte^-^ / d 3 X (\s°^(x),f^(y)\) Suj pa (x) 



(38) 



In order to obtain transport coefficients, a suitable perturbation must be chosen which can be eventually taken out 
from the integral. Physically, it corresponds to the enforcement of a particular hydrodynamical motion and observe 
the response in the stress-energy tensor to infer the dissipative coefficient. The perturbation 5(3 = 1/TSu is taken to 
be a stationary one and non-vanishing only within a finite region V, at whose boundary it goes to zero in a continuous 
and derivable fashion. The perturbation Suj is also taken to be stationary and it can be chosen either to vanish or like 
in Eq. (|2"Tj) ; in both cases, one gets to the same final result. 

Let us then set Slu = and expand the perturbation 5/3 = (0, 0, 5(3 2 (x 1 ), 0) dependent on x 1 in a Fourier series 
(it vanishes at some large, yet finite boundary). Since we want the higher order gradients of the perturbation to 
be negligibly small (the so-called hydrodynamic limit), the Fourier components with short wavelengths must be 
correspondingly suppressed. The component with the longest wavelength will then be much larger than any other 
and, therefore, 5/3 2 can be approximately written, at least far from the boundary, as Asia.(irx /L) where L is the size 
of the region V in the x 1 direction and A is a constant. The derivative of this perturbation reads: 

= yAcos(7rxVi) = ai%(0)cos(7ra; 1 / J L) = <9i#3 2 (0) cos(fca; 1 ) 



where k = n/L. Therefore, by defining k = (A:, 0,0) and plugging the last equation in Eq. ([38 



T 



5T^(y) = lim -&# 2 (0 



e^0 I 



f 1 _ e e(t-t') 

dt 



d x cosk • x{ 



= lim T 9i%(0) Im 



dt 



=.e(t-t') 



T^(x),Tit v (y) 



d 3 xe ,k.x ( T 12 ix))T ^ {y) )q 



(39) 



taking into account that the commutator is purely imaginary. To extract shear viscosity we have to evaluate the 
stress-energy tensor in y = to make it proportional to the derivative of the four-temperature field in the same point 
and we have to take the limit L — > oo which implies V — > oo and k — >• at the same time: 



8T$"(U,0) = lim lim T 9i<58 2 (0) Im 



dt 



1 



Q e(i-t') 



(40) 



where it has been assumed that the integration domain goes to its thermodynamic limit independently of the integrand. 
Because of the time-translation symmetry of the equilibrium density operator po, the mean value in the integral only 
depends on the time difference t—t v . Thus, choosing the arbitrary time t' = t y and redefining the integration variables, 
the Eq. (l4Tfl) can be rewritten as: 



ST^(t v ,0) = lim lim T d^MO) Im 



dt 



dW k ' x < T 1 *(x),Tg v (Q) ) 



(41) 



which shows that the mean value 5Tg V '(t y , 0) is indeed independent of t y , which is expected as 5(3 is stationary. 

We can now take advantage of the well known Curie symmetry "principle" which states that tensors belonging 
to some irreducible representation of the rotation group will only respond to perturbations belonging to the same 
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representation and with the same components 2 . In our case the Curie principle implies that only the same component 
of the symmetric part of the stress-energy tensor, i.e. Ti , will give a non- vanishing value: 



ST^{t y ,0) = lim lim Tdi$3 2 (0) Im / dt 



dW x < T?(x),T?(0) ) 



(42) 



From the above expression, a Kubo formula for shear viscosity can be extracted setting <5/3 = (l/T)Su: 



1) = lim lim Im 



dt 



d^x 



■ik-x / 



Ti 2 (x),T^(0) } 



12/ 



(43) 



which, after a little algebra, can be shown to be the same expression obtained in ref. [9j. Because of the rotational 
invariance of the equilibrium density operator, shear viscosity is independent of the particular couple (1, 2) of chosen 
indices. It is worth pointing out that, had we started from Eq. (|7Tj) instead of Eq. (|70l) . choosing 5ui = or like in 
Eq. (|21[) , we would have come to the same formula for shear viscosity; in the latter case, the third contributing term 
in Eq. (fTTj) would have been of higher order in derivatives of 8(3, hence negligible. 

Now, the question we want to answer is whether equation (|43p is invariant by a pseudo-gauge transformation ([1}, 
which turns the symmetric part of the stress-energy tensor into: 



r; 



(44) 



where: 



(45) 



S being symmetric in the last two indicess. We will study the effect of the transformation on the mean value of the 
stress-energy tensor in the point y = starting from the formula Eq. ([7T|) instead of Eq. ([70]) with 5ui — or like in 
Eq. ([21]) . which allows us to retain only the first contributing term to STg 2 (0). The perturbation 5(3 is taken to be 
stationary and t' is set to be equal to t y = 0. Eventually, the appropriate limits will be calculated to get the new 
shear viscosity. Thus: 



,5T s 12 (0) = 5T] 2 (0) + lim 



dt 



d d x 



d a E al2 (x),d^ 12 (0) ) (d^ (x) +828^ (x)) 



(46) 



lim 



dt 



1-e 



et 



d 3 x < d a E al2 (x),T^ 2 (Q) ) + ( Tj 2 (t,x),d Q S Ql2 (0) >o) + 



We can simplify the above formula by noting that the mean value of two operators at equilibrium can oly depend on 
the difference of the coordinates, so: 



Oi(y),^0 2 (z) ) 



hence, the Eq. (|46| can be rewritten as: 



d 
dx^ 



Oi,0 2 



)o(y - x) 



d 



)o(y - x), 



5T s 12 (0)=5T 1 s 2 (0)-]im 



dt 



d J x 



d 2 



dx a dx' 3 



( E a ' 2 (x),^ 12 (0) )o(9 1 %(x)+9 2 5 / 3i(x)) 



lim 



dt 



d d X 



dx c 



({[^"(^.^(O^o - ([fi 2 (x),S« 12 (0)]) ) (fc%(x) + &# 1 (x)) (47) 



We are now going to inspect the two terms on the right-hand side of the above equation. If the hamiltonian is 
time-reversal invariant, it can be shown (see Appendix B): 



( T|(i,x),S^'(0,O) >o = (-ir°( E ai '(0,0), x) )o = (-l)"°( S^(i, -x), T|(0, 0) ) c 



2 This is true provided that the right hand side of Eq. (1411) is a continuous function of k for k = or that its limit for k — > exists, i.e. 
it is independent of the direction of k 
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where no is the total number of time indices among those in the above expression. Similarly, if the hamiltonian is 
parity invariant, then: 

([i Q «(t,-x),f|(0,0)]) = (-ir<[§^"(i,x),T|(0,O)]) o 

where n s is the total number of space indices. Using the last two equations to work out the last term of Eq. (|47l) one 
gets: 



d 2 



d 3 K 



-2 lim 



dt / d 3 x — 

e J v dx° 



S Ql2 (i,x),^ 12 (0,0) )„(0i%(x)+&#i(x)) 
? 12 (i,x),f' 2 (0,0)l}„(9 1 ft(x)+9 2 ip 1 (x)) (48) 



Now, the two terms on the right hand side of (|48|) can be worked out separately. Using invariance by time-reversal 
and parity, one has: 

([s^'(t,x),S^(0,O)]) o = (-l)"°([s^'(0,O),S m -''(-t,x)]) o 

= (-l)"°([s^'(i,-x),S^(0,O)])o = (-ir + "'([^(t,x),^(0,0)])o = ([s^'(t,x),H^(0,0)])o (49) 

being uq + n s = 6. Hence, the first term on the right hand side of (|48|) can be decomposed as: 



lim 



dt- 



l-e e 



d 3 x 



dt 2 ' 



E 0i ^x), S°«(Q)] ) + 2 J^([H fe « §°« (0) ),, 



+^^T([s feiJ W,S iiJ (0) „,)(W„W,(x)- ( V>;,(x)) 



(50) 



and, similarly, the second term as: 



dt 



1-e 



5 



(9 



d 3 x ^ ( [~ 012 (a,), Tj 2 (0)j )o + ^ < [5 fe12 (*, x) , Tg 2 (0) j ) {d x Sp 2 (x) + (x)) (51) 



- 2 lim 

£-►0 

All terms in Eqs. (|50|) and (|51j) with a space derivative do not yield any contribution to first-order transport coefficients. 
This can be shown by, firstly, integrating by parts and generating two terms, one of which is a total derivative and 
the second involves the second derivative of the perturbation <5/3. The total derivative term can be transformed into 
a surface integral on the boundary of V which vanishes because therein the perturbation 5/3 is supposed to vanish 
along with its first-order derivatives. The second term, involving higher order derivatives, does not give contribution 
to transport coefficients at first order in the derivative expansion. Altogether, the Eq. (|48[) turns into: 



5T^(0) = 5T^(0)-lim 



dt 



1 -e e 



d 3 x5. 



-2 lim 

£->0 



dt 



1-e 



et 



■IV 



2 ui2 (z),S ui2 (0) >o(&%(x) +3 2 $i(x)) 



d 3 x5 t ( E ul2 (x),T s 2 (0) ) (d 1 Sf3 2 (x) + d 2 5(3 1 (x)) + O(d 2 6f3) 



which can be further integrated by parts in the time t, yielding: 



(52) 



ST' s 12 (0) = 6T^ 2 (0) - lim / dt (S(t) - £C £t ) / d 3 x ( S ul2 (x), S ul2 (0) ) (5 1 ^ 2 (x) + &#i(x)) 



-2 lim 



dt c e 



d 3 x 



E 012 (x),T s 2 (0) ) (9i%(x) +9 2 %(x)) +0(d 2 S(3) 



(53) 



provided that, for general space-time dependent operators 0\ and C* 2 

d 



lim 



d 3 x e n£t 



dt 



( Oi(i,x),O 2 (0,0) ) =0 



3,_ n net i 



lim / d x e 

t— > — oo 



Oi(i,x),O 2 (0,0) ) = 
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with n = 0,1, which is reasonable because thermodynamical correlations are expected to vanish exponentially as a 
function of time for fixed points in space 3 . 

From Eq. (|53p the variation of the shear viscosity can be inferred with the very same reasoning that led us to 
formula (T4"3"l). that is: 



An = T) — r\ 



lim lim Im 

£->0 fe^O 



2 lim lim Im 

e->0 fc->u 



dt (5(t)-ee et ) 



d 3 x e lkxl i 



~^(t,x),~ u -(0,0) ) 









f d 3 x e ikxl ( 


/ —oo J 





If the first integral is regular, then the e — > limit kills one term and the (f54"f reduces to: 

d 3 x cos kx 1 ( [S 012 (0, x), S 012 (0, 0)1 ) 



An = rj — n = — lim 

fe^O 



2 lim lim Im 

E->0 fc^O 







/ dt e et 


f d 3 x e ikx± ( 


t — oo J 





(54) 



(55) 



In general, this difference is non- vanishing, leading to the conclusion that the specific form of the stress-energy tensor 
and, possibly, the existence of a spin tensor in the underlying quantum field theory affects the value of transport 
coefficients. The relative difference of those values depends on the particular transformation ([1]), hence on the 
particular stress-energy tensor. In the next Section a specific instance will be presented and discussed. 

An important point to make is that the found dependence of the transport coefficients on the particular set of 
stress-energy and spin tensor of the theory is indeed physically meaningful. This means that the variation of some 
coefficient is not compensated by a corresponding variation of another coefficient so as to eventually leave measurable 
quantities unchanged. This has been implicitely proved in Sect. IIVI where it was shown that total entropy itself 
undergoes a variation under a transformation of the stress-energy and spin tensor (see Eq. (|3"5j0 . 



VI. DISCUSSION AND CONCLUSIONS 



As a first point, we would like to emphasize that in our arguments spacetime curvature and gravitational coupling 
have been disregarded. On one hand, this shows that the nature of stress-energy tensor and, possibly, the existence of 
a fundamental spin tensor could, at least in principle, be demonstrated independently of gravity. On the other hand, 
for each stress-energy tensor created with the transformation ([T]), it should be shown that an extension of general 
relativity exists having it as a source, which could not be always possible. 

An important question is whether a concrete physical system indeed exists for which the transformation ([T]) leads to 
actually different values for e.g. transport coefficients, entropy production rate or other quantities in nonequilibrium 
situations. For this purpose, we discuss a specific instance regarding spinor electrodynamics. Starting from the 
symmetrized gauge-invariant Belinfante tensor of the coupled Dirac and electromagnetic fields, with associated S — 0: 

= l - (% 7 ^ v'^ + V^*) + F\F Xv + ^ U F 2 (56) 

where = — ieA^ is the gauge covariant derivative, one can generate other stress-energy tensors with suitable 
rank three tensors and then setting $ = —S' where S 1 is the new spin tensor, according to ([1]). One of the best known 
is the canonical Dirac spin tensor: 

({ } stands for anticommutator) which is gauge-invariant and transforms the Belinfante tensor (|56")l back to the 
canonical one obtained from the spinor electrodynamics lagrangian (see also [l[ for a detailed discussion). However, 
this is totally antysimmetric in the three indices A, //, v and thus the variation of T operator (see Eq. (|2"2"]) as well as 



3 There might be singularities on the light cone, however for fixed x and and integration over a finite region V, in the limit t — * — oo 
light cone is not involved 
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transport coefficients, which depend on the symmetrized S tensor (|45|) vanish. Nevertheless, other gauge-invariant 
<&-like tensors can be found. For instance, one could employ a superpotential: 

8m V / 8m V / 

which is the gauge-invariant version of the one used in ref. [Toj ] to obtain a conserved spin current. This superpotential 
gives rise to a non- vanishing spin tensor as well as a S tensor (see Eg. 

lorn V / 

hence a variation of thermodynamics. By noting that the structure of the above tensor is very similar to the Belinfantc 
stress-energy tensor (|56|) . it is not difficult to find a rough estimate of the variation of e.g. shear viscosity induced 
by the transformation. Looking at Eq. (|55p we note that S 012 mainly differs from T 012 in Eq. ([5B| by the factor 
1/m. The last term on the right hand side of Eq. (|5S|) tells us that the dimension ofH is that of a stress-energy tensor 
multiplied by a time, and therefore this term must be of the order of rjh/mc 2 T where r is the microscopic correlation 
time scale of the original stress-energy tensor or the collisional time scale in the kinetic language and r\ the shear 
viscosity obtained from the original stress-energy tensor. Thus, the expected relative variation of shear viscosity from 
Eq. ([55)1 in this case is of the order: 



V 

which is (as it could have been expected) a quantum relativistic correction governed by the ratio (A c /c)/t, A c being 
the Compton wavelength. For the electron, the ratio A c /c « 1CP 21 sec, which is a very small time scale compared to 
the usual kinetic time scales, yet it could be detectable for particular systems with very low shear viscosity. 

It is also interesting to note that the "improved" stress-energy tensor by Callan, Coleman and Jackiw [l3] with 
renormalizable matrix elements at all orders of perturbation theory, is obtained from the Belinfante's symmetrized 
one in Eq. ([55| with a transformation of the kind ([1]) setting (for the Dirac field and vanishing constants [Oj]): 

6 

and requiring §' = § = so that $ A '^ = d a Z aX ^, hence: 

$a,^ = __ (g^d* 1 - g^d") 
6 

GV" _ ifS^iAf i S^^^m = —Z 
2 V ' 6 

6 

which is just the improved stress-energy tensor 13]. It is likely (to be verified though) that the aforementioned 
modified stress-energy tensors imply a different thermodynamics with respect to the original Belinfante symmetrized 
tensor. 

To summarize, we have concluded that different quantum stress-energy tensors imply different values of nonequilib- 
rium thermodynamical quantities like transport coefficients and entropy production rate. This reinforces our previous 
similar conclusion concerning differences of momentum and angular momentum densities in rotational equilibrium [lj] . 
The existence of a fundamental spin tensor has, thus, an impact on the microscopic number of degrees of freedom 
and on how quickly macroscopic information is converted into microscopic. The difference of transport coefficients 
depends on the particular form of the tensors and in the examined case it scales like a quantum relativistic effect 
with h/c. Therefore, at least in principle, it is possible to disprove a supposed stress-energy tensor with a suitably 
designed thermodynamical experiment. 



2 



14 



Acknowledgments 

We are grateful to F. Bigazzi, F. W. Hehl and D. Seminara for useful discussions and suggestions. 



References 



[1] F. Becattini and L. Tinti, Phys. Rev. D 84, 025013 (2011) 

[2] F. Halbwachs, Theorie relativiste des fluids a spin, Gauthier-Villars, Paris, (1960). 
[3] F. W. Hehl, Rept. Math. Phys. 9, 55 (1976). 

[4] D. N. Zubarev, Sov. Phys. Doklady 10, 850 (1966); D. N. Zubarev and M. V. Tokarchuk, Theor. Math. Phys. 88, 876 

(1992) [Teor. Mat. Fiz. 88N2, 286 (1991)]. 
[5] H. A. Weldon, Phys. Rev. D 26, 1394 (1982). 
[6] F. Becattini, Phys. Rev. Lett. 108, 244502 (2012) 

[7] D. N. Zubarev Non Equilibrium Statistical Thermodynamics, Nauka, Moscow (1971). (english translation: New York, Con- 
sultant Bureau, 1974); D, N. Zubarev, V. G. Morozov, G. Ropke, Statistical Mechanics of Nonequilibrium Processes. 
Volume 2: Relaxation and Hydrodynamic Processes, Akademie Verlag, Berlin (1997). 

[8] R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II. Nonequilibrium Statistical Mechanics, Springer- Verlag, Berlin 
(1985). 

[9] A. Hosoya, M. Sakagami and M. Takao, Annals Phys. 154, 229 (1984). 
[10] S. R. De Groot, W. A. van Leeuwen, Ch. G. van Weert, Relativistic kinetic theory, North Holland (1980). 
[11] R. Balian, From microphysics to macrophysics, Springer, Heidelberg (2007). 
[12] L. Landau, L. Lifshitz, Statistical Physics, Pergamon Press (1980). 
[13] C. G. Callan, Jr., S. R. Coleman and R. Jackiw, Annals Phys. 59, 42 (1970). 



APPENDIX A - Relativistic linear response theory with spin tensor 



We extend the relativistic linear response theory in the Zubarev's approach to the case of a non-vanishing spin 
tensor. The (stationary) nonequilibrium density operator is written in Eq. (|13[) . with T expanded as in Eq. (| 14[) . 
As has been shown in Sect. UH at equilibrium, only the first term of the T operator survives in Eq. (|14[) ; therefore, 
one can rewrite that equation using the perturbations 5/3, St; and Soj which are denned as the difference between the 
actual value and their value at thermodynamical equilibrium: 



dt e e(t -*' } J dSn, (f^Sp^x) - ?<J£(a:) - \§ 1 '^ 5oo^{x) 



lim 

e^0 



..... / J 

2 £->o 



1 lim / dt 



x e^ 4 "*'' (f^id^ix) + d^ v {x)) + f^(d^ u (x) - d M 50„(x) + 2&v(a:)) 

-S^dxSuj^ix) - 2fd^{x)) (57) 
where it is understood that x = (t, x). 

In fact, we will use a rearrangement of the right-hand-side expression which is more convenient if one wants to work 
with an unspecified, yet small, Suj. Therefore, the above equation is rewritten as: 

1 

T 

lim J' dt e^*" ^ J d 3 x (t 0v 5I3 v (x) - ±S°^6u u „(x) -?<5^)) (58) 



d 3 x ( T°^„(t',x) - ?£(i',x) - -S^avCi'.x) 



what it can be easily obtained from Eq. (|13j) integrating by parts in time. 

For the sake of simplicity we calculate the linear response with £ cq = <5£ = 0, but it can be shown that our final 
expressions hold for £ cq ^ (in other words with a non- vanishing chemical potential fi ^ 0). Let us now define: 



A=- Jd 3 x (T ^„(i',x) - i<S°^„(i',x) 
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and: 



so that: 



B = hmj dt e E(t - 4,) ^ J d 3 x (f 0v 6f3„(x) - l^Su^a 



P=|exp[-T] = icxpIA + S] 



(59) 



with Z = tr(exp[A + B])s. 

The operator B is the small term in which p is to be expanded, according to the linear response theory. It can can 
be rewritten in a way which will be useful later on. Since: 



d'x -|-/° = /,i-\«v(.f ""(,•), <>•„(•/•.)) - / d : x ( v,T"',.,.) a ;,.,(,■! 



Jay 



then: 



B = -faaJ dte'W J d 3 x {t^SP^x) - ~ (^'""^(a;))) - J^dS n z f™(x) 5/3 u (x) 
The perturbation <5/3 must be chosen such that (5/3 1 ay = so that only the bulk term survives in the above equation: 



At the lowest order in B: 

Z = ti{e A+S ) ~ tr(e^ 
and, according to Kubo identity: 



1 + B 



) = z LE (i + (S) LE ) => i ~ J-(i - (s) LE ) 

^ ^LE 



1+ [ dze< US )Be- zA 
Jo 











1 + 











dz e ZA Be 



(60) 



(61) 



(62) 



where the subscript LE stands for Local Equilibrium and implies the calculation of mean values with the local 
equilibrium density operator (see Sect. IIV[) . Thereby, putting together (I5T|) and (p32"|) and retaining only first-order 
terms in B: 



1 - (B)le ) Ple + I dz e zA Be- zA p hE , 



hence the mean value of an operator 0(y) becomes: 



(0(y)> ^ (l - (-B)leJ (0(2/))le + (O(y) jf dz e^Be"' 

Let us focus on the last term, which, by virtue of (IBUl) . contains expressions of this sort: 

(d(x)X'(z,t,x)) LE = (6(x)e zX X(t^)e~ zX ) LE 
where X stands for components of either T or S or d$S. From the identity: 

(d{y)X'(z,t,x)) Jj - E = [ dr (6( 2 /)d r A > '(z,T,x)) LE + lim (d(y)X'(z, r, x)) LE , 



(63) 
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and the observation that correlations vanish for very distant times (check footnote [3]) , one obtains: 

(d(y)X'(z,t,x)) L E= / dr (d(y)d T X'(z,T,x)) LE + lim ((5( 2 /))le(X(t, x)) LE , (64) 



where we have also taken advantage of the commutation between expL4] and exp[±z/l]. 

We now approximate Q the local equilibrium density operator with the nearest equilibrium operator p in Eq. (|26p , 
which also implies that: 

A ~ —H/T 

where H is the hamiltonian operator (which ought to exists given the chosen boundary conditions). The straightfor- 
ward consequence of this approximation is that the second term on the right hand side in Eq. (|64[) can be written 
as: 

(X(-oo,x)) LE ~ (X(-oo,x)) = (X(t,x)) 

because the mean value is stationary under the equilibrium distribution. Therefore, the Eq. can be approximated 
as: 

(d(y)X'(z, t, x))lb * I dr (d(y)d T X'(z, r, x)) + (6(y)) Q (X(t, x)) , (65) 

J — CO 

and the as: 

(d(y)) ~ (i_(B) )(6(y))o+ / d* (O(y) e-^/ T B e zl? / T }o (66) 



Once integrated, the second term in (|65|) gives rise to a term which cancels out exactly the (B}o(0(y))o in the equation 
above, which then becomes: 



(0(y)) ~ (O(y)) + dz dr {0(y)d T X'(z, r, x)) (67) 

JO J-oo 

Let us now integrate the last term on the right hand side in z: 

[ dz [ dr(6(y)a T X'(z,r,x)) = i [ du [ dr (6(y)d T e- uS X(r, x)c u % 

JO J-oo P JO J-oo 

where ft = 1/T and (3z = u. As H is the generator of time translations: 

i / du [ dr {d(y)d T e- u "X( T , x)e u % = \ [ du [ dr (d(y)d T X(r + iu, x)) 

P JO J-oo P JO J-oo 

= A / dr du ^ {(O(y)X(r + iu, x)) ) =^J_ ((O(y)X(r + 0, x)) - (d(y)X(r, x)) 
On the other hand: 

(O{y)X(T + iP,x)) = tr(p O(y)e-^X(r,x)e+^) = i-tr( e -^0(y)e-^X(r,x)e^) 

^o 

= i-tr(0(y)e-^X(r,x)) = tr(X(r, X ))p d(y)) = (X(r, x)O(y)> 
^o 



Hence, putting the last three equations together, we have: 



1 '« 



dz / dr(O(y)a T X'(z,r,x)) = ^ / dr ([X(r,x), O(y)]) (68) 
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Substituting now X with its specific operators, Eq. (|67|) can be expanded as: 



S(d(y)) = (d(y)) - (O{y)) * lim j* I dt / dr / d 3 x ( f ""(r, x), O(y) ^^(a) 



ilimi dte^*"*'^ dr / d 3 x 



: lim i / dt e e (*-*'> / dr / d 3 x 

£->0 



-— lim 

-— lim 

2 e ->o i/3 



— oo 

t' 



r^(r,x),0(y) ) 3 M <S3„(a;) 



dt e 



s(t-t') 



d 3 x 



5 ''"'(t,x),O(i/) >o&v0»0 



die^*"*') / dr / d 3 x<[<^>,x),6>)])ojU^(*) 



The first term on the right hand side of the above equation can be integrated by parts using: 



dte £ (*-*') / dr/(r 



dt — 
at 



d /e £ (*-*')\ I* 



dr /(r) 



dr /(r) 



t' e e(t-t') 
dt - 



/(f) = / dt 



t' l _ e s(t-t') 



fit) 



so that the Eq. (|69|) can be finally written: 



S(0(y)) = lim - 
lim 

2 e^O 

lim 

' e^o i/3 



1 /•*' l-e^*-*') 



dt 



d 3 x([T^(x),O(y)j) d^(x) 
'5 D ^(x),O(y)l) <5^(x) 



die^"*') / dr f d 3 x([5 n ^(r,x),6(y)]) ^^(x) 



dt e £ ('-*') / d 3 x 



(69) 



(70) 



Another useful (equivalent) expression of 6(0(y)) can be obtained starting from the expression (ITi)) of T, where 
the continuity equation for angular momentum is used from the beginning. Repeating the same reasoning as above, 
it can be shown that one gets: 



5(0(y)) = lim -5 



t' 1 _ e e(t-t') 

dt 



d 3 x 



lim 

e^o 2i(3 



1 ft' l- e e(t-t') 



dt 



d 3 x 



T^(x),0(y)\ )o(dMx) + duSPtix)) 
Tf (z), 6(y)] ) (dMx) - duSP^x) + 2Su lu ,(x)) 



1 l im L / dt e^*"^ / dr / d 3 x ( S x ^{r, x), 6(y) ) d x 8w^{x) 



2 e^O if] 

As we have pointed out, these expressions hold when po has a non-vanishing chemical potential. 



(71) 



APPENDIX B - Commutators and discrete symmetries 

We want to study the effect of space inversion and time reversal on the mean value of commutators like: 

([5f •"»(i J x),0^-"-(0,0)]>o 

where 0% and O2 are tensors of rank respectively m and n. 

The equilibrium density operator p = cxp[—H/T]/Z is symmetric for space-time translations and rotations, as 
well as time reversal and parity if the hamiltonian is itself parity and time reversal invariant. The symmetry under 
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this class of transformations allows to simplify the above expression. For any linear unitary transformation U which 
commutes with p one has: 

(6)0 = tr (po 6) = tr (iHpbU 6) - tr (p DOIT 1 ) = (UOU-^o 
Taking U = T(a) with T(a) a general translation operator: 

( [d^ """ (t, x), O?-"" (0, 0)] ) = < [6J* 1 ' """ (t + a°, x + a), 6?-"» (a , a)] ) 
and so, setting (a , a) = (— t, — x): 

( (t, x), d^ - Vn (0, 0)] ) = ( [6r (0, 0), O^-"" (-t, -x)] } 
Similarly, for a space inversion: 

( [Of (i, x), O^ 1 (0, 0)1 } = (-l) n =+ m = ( fe 1 '"""* (t, -x), Of" (0, 0)1 )o 

where m s and n s are the number of space indices among fii, - • • p m and v\, ■ ■ ■ v n respectively. 

The time reversal operator 6 is antiunitary, thus a point-dependent scalar operator A(t, x) transforms as follows: 



el^x)^- 1 = lt(-t,x) 



whence, for commutators: 



9 



9" 1 



Bt(_ t ,x),j4t(-t,x) 



For hermitian operators what is then changed is the order of the operators beside their time argument. For tensor 
hermitian observables and time-reversal symmetric hamiltonian, one obtains: 

( [Of (t, x), d?-"» (0, 0)] ) = (-l)™o+no ( (0, 0), Of 1 -"™ (-t, x)] ) 

where m and n are the number of time indices among \x\ , • • • /z m and v\ , ■ ■ ■ u n respectively. 



